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DYNAMIC ANALYSIS OF NONCONTACTING FACE SEALS 

by I, Etsion 

ABSTRACT 

The dynamic behavior of a noncontacting’ coaed face seal is analyzed 
taking into account various design parameters and operating conditions. The 
primary seal ring motion is expressed by a set of nonlinear equations for 
three degrees of freedom. These equations, which are solved numerically, 
allow identification of two dimensionless groups of parameters that affect 
the seal dynamic behavior. Stability maps for various seals are presented. 
These maps contain a stable-to-unstable transition region in which the ring- 
wobbles at half the shaft frequency. The effect of various parameters on 
seal stability is discussed and an empirical expression for critical stability 
is offered. 
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NOMENCLATURE 
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(1 - R)/H^H''(1 - R.) 
linearity constant, eq, (38) 
seal balance ratio, (r^ - r^)/(i*Q - r?) 
seal center- line clearance 
equilibrium center-line clearance 
axial force 
closing force 
initial spring load 
dimensionless force, F/Sr^ 
dimensionless film thickness, h/C 
dimensionless minimum film thickness 
film thickness 

ring mass moment of inertia about a diameter 

2 4 

dimensionless moment of inertia, I«o.^cySr’ 
spring constant 

dimensionless spring constant, K*/Sr^ 
moment 

external moment 
dimensionless moment, M/Sr' 
ring mass 

2 2 

diriensionless mass, m* 0 J^C^/Sr^ 
number of seal revolutions 
dimensionless pressure, p/s 
pressure 
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t* 
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a* 

a 
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Subscripts : 


dimensionless radius, r/r 

radial coordinate 

seal balance radius 

ring radius of gyration 

seal parameter, 6Mco(r^/C^) (1 - R.) 

time 

dimensionless time, cot* 
axial displacement 
dimensionless displacement, Z*/C^ 
tilt angle 

normalized tilt, o!*r_/C 
coning angle 

normalized coning, /3*r^/C^ 
nutation 

normalized nutation, 
coning parameter, jS*r^/C 
tilt parameter, y*r^/C 
angular coordinate 
viscosity 
precession 

shaft angular velocity 


1,2,3 

d 

i 

m 


axes 1, 2, or 3, respectively 
hydrodynamic 
inner radius 
mid radius 


o 


outer radius 
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s hydrostatic 

sp springs 

X, y, z axes, x, y, or z, respectively 

INTRODUCTION 

The presence of a fluid film between the mating faces of the primary seal 
ring and seal seat in a mechanical face seal (fig, 1) was established in the 
early 1960’s by Denny [1]. Since then several different mechanisms which 
provide sources of fluid film load support have been proposed. These mech- 
anisms include surface waviness, angular misalinement, thermal and mech- 
anical distortions, and phase change. The state of the art is well summa- 
rized in reference [2] and some of the latest results on seal theory are 
presented in references [3-9]. From all these works it is evident that the 
pressure developed in the lubricating film is strongly dependent on the shape 
and magnitude of the film thickness between the mating faces. Almost all the 
analyses published so far assume a given fixed shape of the film thickness 
and then calculate the resulting pressure distribution from which forces, 
moments, flow rates, and friction can be derived. In reality, however, the 
elements forming the seal are not rigidly mounted. Hence, the film thick- 
ness is both space and time dependent. 

The main objective in any seal is to minimize leakage. To this end the 
separation of the mating faces should be as small as possible. However, too 
small a clearance increases friction losses and can result in excessive wear. 
Therefore a proper sealing gap has to be maintained at all times to ensure 
safe and reliable operation. Short life, due to surface damage, and excessive 
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leakage are common problems in face seals, These problems can be at- 
tributed to unstable operation which causes variations of the sealing gap be- 
yond permissible limits. Indeed, dynamic instability in the form of vibra- 
tion of flexibly mounted elements was experimentally observed in face type 
seals [10 - 14]. 

The primary seal ring in a face seal (fig. 1) is flexibly mounted so 
that only its circumferential rotation is prevented. Hence, the ring has 
five degrees of freedom; axial displacement, two angular tilts about two 
orthc^onal diameters, and two perpendicular radial displacements. Such 
a system can be vulnerable to several sources of instabilities which are 
discussed in the seal literature mentioned above. In order to avoid unstable 
operation a mechanical face seal has to be designed with a great care re- 
garding the mass, inertia, stiffness, and damping combinations for selected 
geometries and operation conditions. To this end a complete dynamic 
analysis is required which takes into account the various parameters affect- 
ing seal stability. 

Only a few attempts have been made so far to analyze seal dynam- 
ics. Most of these analyses deal with only one degree of freedom or over- 
look some or all of the fluid film effects [15 - 19]. A proper treatment of 
the problem is presented in reference [20] but it treats a special case 
where the system forces and, moments are provided by a supporting gas 
bearing and not by the sealir^ dam. In more recent work [21, 22] the 
motion of a flexibly mounted ring in a noncontacting face seal is described 
in its three major degrees of freedom (one axial and two angular). Hydro- 
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dynamic, hydrostatic, and squeeze film effects as well as the contribution 
of the flexible support to the system dynamics are considered and it is 
shown how the seal stability is affected by various design parameters. The 
analysis in references [21] and [22] is, however, somewhat limited in that 
it is based on small perturbation in the seal clearance and assumes flat 
faces. 

The objective of this paper is to present a more comprehensive dynamic 
analysis of a noncontacting face seal. The nonlinear equations of motion of 
the primary seal ring in its major three degrees of freedom are solved 
numerically, giving the complete transient behavior of the ring. In addi- 
tion, coning of the mating faces is included, thus, resulting in a more real- 
istic model in which balance ratio and initial spring loading are also con- 
sidered. 

THEORY 

The theoretical model is shown in figure 2. The seal seat is parallel 
to the plane xy of an inertial reference xyz and is rotating at a constant 
ai^ular velocity o) about the z axis. The primary seal ring has three deg- 
rees of freedom, it can move axially along the z axis and tilt about the x 
and y axes. A rotating coordinate system 123 coincides with the principal 
axes of the ring so that axis 3 is perpendicular to the plane of the ring. 

The coordinate system 123 rotates in the inertial reference xyz so that 
axis 1 always remains in the plane xy and axis 2 is directed to the in- 
stantaneous point of maximum film thickness. Thus, the orientation of 
coordinate system 123 in the inertial reference xyz can be defined by the 
nutation angle y* measured from axis z to axis 3, and the precession 
angle i// measured from axis x to axis 1. 


7 


Due to the tiny clearance C (of the order of few micrometers in typical 
seals) the nutation angle y* is also very small. Hence, as shown in [21] 
the equations of motion of the seal ring are 


F, = m Z 




( 1 ) 

( 2 ) 

(3) 


where a* and a* represent the acceleration of the ring about the x and 

J 

y axes, respectively. 

The system of force F* and moments M* and M* is composed of 
fluid film pressures, pressure from mechanical springs, and sealed pressure 
on the back of the ring. The moments acting on the ring are more readily 
obtained with respect to axes 1 and 2 of the rotating coordinate system 123 
and then transformed into the inertial reference xyz. This transformation 
takes the form of (see fig. 3) 


M* -- Mj cos - Mg sin \}^ (4) 

My = Mj sin \(y + Mg cos ( 5 ) 

The moment Mg contributes the terms Mg sin y*sin and 

-Mg sin y * cos to M* and M*, respectively. However, using Euler's 

equations of motion it can be shown that for y * « 1 , which is our case. 

Mg itself is of order y*2 and therefore can be neglected. 
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The moments 



Mj" and Mg are calculated from 
pr^ cos 0 dr d0 - M* 


( 6 ) 




sin Q dr d0 


and the force F* is 

iZi 



( 7 ) 


( 8 ) 


where M* is the restoring moment provided by the flexible support (e.g. , 
mechanical springs) and F* is a closing force contributed by both the 
flexible support and the pressure balance of the seal (see appendix 1). 

The pressure distribution, p, is found from a solution of the Reynolds 
equation which, by usii^ the narrow seal approximation [23], reduces to 

J./h^ = 6/u (cu — + 2 (9) 

3r\ ar/ \ 30 3t*/ 


The boundary conditions of equation (9) are 


p = p. at r = r. 
p = Po at r = Tq 

Because of its linear nature, equation (9) can be solved separately for the 
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hydrostatic, hydrodynamic, and squeeze components of the pressure 
[6, 7, 8], The hydrostatic pressure component is 


Ps = Po 


(Pq - Pi) 




( 10 ) 


The hydrodynamic and squeeze pressure components can be combined to give 


Pj = -3 — + 2 

“ '30 at=*-y 


(r^ - r)(r - r.) 


m 


( 11 ) 


and the total pressure is the sum of Pg and p^. Hence, the complete solu- 
tion of equation (9) is 

P = Ps + Pd (12) 

Assuming that the mating faces have some coning which can result 
from either thermal and mechanical distortion, or can be intentionally ma- 
chined into one of the surfaces, the film thickness distribution h is (see 
fig. 2) 

h = C + y*r cos 6 + ^*(r - r.) (13) 

Visualizing an angle cp measured from the y axis to a plane defined by the 
z axis and a point (r, 9) on the rir^ face (fig. 2), we have for y* «1 

e = cp - \{/ ( 14 ) 

where <p is fixed in time. Hence, substituting (14) in (13), differentiating 


and neglecting curvature effects we have 


— = -y * r^ sin 0 (15) 

de ^ 

and 

JJL = z* + y * r cos 9 + sin 9 (16) 

3t* m ^ m 

where the dots indicate differentiation with respect to time. 

Using the normalized parameters y = y* R = r/r^, t = wt*, and 

Z = zVCq where is the designed clearance (see appendix 2). Noting 
that 

C = + Z* = C^(l + Z) 

and normalizing the film thickness in the form 

E = ^= 1+ €R cos 0 -!- 6 (R - R.) 

C ^ 

we have, from equations (10) and (11) , 

Ps = Po - <Po - ^ H) I 


(17) 


(18) 



r sin 6 ~ Z ~ y cos 0 


A ^ ~ ^i 
(1 + Z)2 ^ " ^i 


(19) 


In equations (18) and (19) A is simply 
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1 - R 


HmH^d - Ri) 

and P is a dimensionless pressure obtained from P = p/S. The term S 
is a seal parameter having -dimension of pressure and defined as 


S = 6 yw (1 - Rj)2 (21) 

\^ol 

Using the seal parameter S to normalize forces and moments in the form 
F = F*/Sr^ and M = M*/Sr^, recalling the normalized time t and axial 
displacement Z, and using a normalized tilt in the form a - the 

equations of motion (1), (2), and (3) become 


F^ = mZ 
z 




My = toy 


Equations (22) to (24) can be solved for the accelerations from which the ring 
motion can be found. 

NUMERICAL CALCULATION 

In order to examine the dynamic behavior of the seal the flexibly 
mounted ring is disturbed from its equilibrium position y = y = 0, Z = Z = 0. 
The new ,tilt and conirig parameters e and 6 are calculated from 
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e = y+!o = ._JL_ (25) 

C 1 + Z 

5 = ^ — (26) 

C 1 + Z 

These parameters are used in equations (17) and (20) to find the film thickness 
distribution H and the parameter A, respectively. The ring area is divided 
radially and circumferentially by a fine grid and the pressures Pg and 
are calculated at each grid point usir^ equations (18) and (19). The total pres- 
sure P is then calculated at each point by adding Pg and P^j, hence, 

P = P„ P , (27) 

s d 

Any negative value of the total pressure P is replaced by P = 0, thereby 
taking account of possible cavitation. The discrete pressures are then inte- 
grated numerically to obtain the fluid film force and moments with the cavi- 
tation effect included. Hence, from equations (6), (7), and (8), 



( 29 ) 



( 30 ) 


The moments and Mg are used to calculate M^^ and M by equations 

• • • • « t 

(4) and (5) and the accelerations Z, a^, and ofy are calculated from equa- 
tions (22), (23), and (24). The new velocities Z, a and a and new dis- 

X y 

placements Z, and ot^ are found by a time integration routine and the 
new nutation y and precession t// are calculated using the relations (see 


fig. 3) 


01 = y cos \}j 

A 

(31) 

Of - r sin 4^ 

(33) 

which are valid for small angle.<3. Differentiating equations 
with respect to time we have 

(31) and (32) 

• * « 

oi = y cos xl' - 74^ sin \0 

(33) 

• f • 

Qfy = y sin 4^ + yif^ cos ;// 

(34) 

• • 

Hence, y and ;// can be calculated from the already known a^, a , y, 

y 

and i//. A new film thickness H and pressure distribution P can now be 
found and the whole process is repeated. 

The procedure described above is carried on in time giving the tran- 
sient dynamics of the seal. At each time step the minimum film thickness 
is calculated and tested against a failure criterion. As can be seen from 
figure 2, contact can take place on either the inner or the outer radius of 
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the risig at 0 ^ tt. Hence, by equation (17), the minimum film thickness is 
calculated from either 


^min = 1 ^ 


(35) 


or 




(36) 


whichever is smaller. The transient behavior of serves as an indica- 
tion of seal stability. Starting from the instant of the disturbance, can 

decrease, increase, or reach a constant value. If increases or reaches 

^min " const the seal is considered stable (with the special case of = 1 


for parallel opei’ation). If, however, is ever decreasing, and even- 

tually diminishes,' a face contact occurs and the seal is unstable. 

RESULTS AND DISCUSSION 

Three different values of seal radius ratio were examined, these 
are: 0.8, 0,9, and 0.98. At each radius ratio a variety of design parame- 
ters and operation conditions were tested and three distinct modes of opera- 
tion were identified. The data are presented in figures 4, 5, and 6 in the 
form of stability maps. An interesting and useful result is the relation found 
between the dimensionless groups of parameters (rg/rgp)^m*co^/K* and 

^^m'^^sp^ ^^o “ Pi^ ^ ^'o^^*^o> ''^here i-g is the ring radius of gyration 
relating mass and inertia by 


It = i m*r^ 

9 o 


From the many cases run on the computer, it was found that the critical sta- 
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bUity, or tlie transition from stable to unstable mode of operation, is char- 
acterized by the expression 



+ 


a 





K*C„ 



(37) 


o 

with a as a linearity constant. Whenever m*(o /K* is less than the 

p 

critical value (m*cu given in equation (37), any disturbance of the 

ring from its equilibrium position vanishes after awhile and the seal is 

p 

stable. When m^covK* is larger than the critical value, the seal becomes 
unstable, any disturbance increases in time causing eventual seal failure 
due to face contact, The transition from stable to unstable operation takes 

p 

place in a narrow band adjacent to the critical value (m*o> /K*) as is in- 
dictated by the shaded areas in figures 4 to 6. In this transition region the 
nutation y* and, hence, the minimum clearance, reach a constant value 
while the ring wobbles at half the shaft frequency. 

The linearity constant, a, in equation (37) is a function of the radius 
ratio Rj and the normalized coning | 3 . With the exception of the cases cor- 
responding to Rj = 0.8 and |3 <5, the constant a in all the other cases was 
found to be approximated fairly by 


9 1 “ I^R. 

a = -STTRj'd - Rj)^ ^ 

2 + |3(1 - R^ 


(38) 


This constant, which represents the slope of the stability threshold lines in 
figures 4 to 6, is maximized, as can be seen from equation (38), whenever 
the coning is 


(39) 


B .= S 

R.(l - R.) 

Hence, at each radius ratio the stability region can be increased by increas- 
ing the normalized coning up to the optimum value of Further increase 

of the normalized coning results in reducing the stability. 

Ten representative cases at R^ = 0.9 were selected to demonstrate the 
various modes of seal behavior and the affect of various parameters on its 
stability. The characteristics of these cases are summarized in table I. 

Some of the computer program results are shown in figures 7 to 11, and the 
effects of various parameters are discussed below. 

p p 

Dimensionless Group (r /r ) m*o) /K* 

o 

' ' •' ' p p 

The effect of values of the group of parameters (r /r ) m*o) /K* on 

O ®lr 

seal stability is shown in figures 7, 8, and 9, Figure 7 presents time varia- 
tion of the minimum clearance, H_. , vs, the number of seat revolutions 
after the ring is disturbed from its equilibrium. Figures 8 and 9 present 
the motion of the rii^ in its angular and axial degrees of freedom, respec- 
tively. Starting with figure 7 and using the stability map for R^ = 0, 9 (fig. 5) 
we see that case 1 is in the unstable regime, case 2 in the transition region, 
and cases 3 and 4 are in the stable regime of operation. From figure 7 we 
see that in case 1, after the ring is slightly disturbed, the minimum clear- 
ance is slowly decreasing and after about 20 revolutions the seal fails due to 
face contact. The failure in case 1 is associated with increasing amplitude 
of angular vibrations as described by 0 !^ in figure 8 and also an increase in 
the center-line clearance, C, as shown by the axial displacement in figure 9. 
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In case 2 the initial disturbance disappears rapidly, after about 1,5 
revolutions. From there on the minimum clearance is maintained constant 
at = 0. 51 as shown in figure 7, while of (fig. 8) oscillates at a con- 

stant amplitude and at half the shaft speed. In the absence of an external 
moment the ring rotation about the x axis is identical to the one about the y 
axis with only a phase shift between the two motions. In the transition be- 
tween stable and unstable modes of operation, represented by case 2, the 
phase shift is 90 degrees, hence, and result in a constant nutation 
y and a constant precession i// = w/2. The axial displacement in case 2 
(not shown) also maintains a constant value which is Z = 0, 145, Thus, the 
center- line clearance is increased slightly due to the constant nutation y. 

Case 3 is a stable case witli a rapid decay of any disturbance as shown in 
figures 7 and 8. Case 4 is similar to case 3 with the addition of an external 
moment M^. As a result of tlie external moment, which can be caused, for 
example, by manufacturing tolerances in the supporting springs, the equili- 
brium position of the ring is not parallel to the seat. This can be seen from 
figure 7 where H^. reaches a constant value of 0. 29 rather than 1.0 as 
in case 3. Tlie ring remains tilted with respect to the seat, with o' and 

A 

Q!y also maintaining constant values (not shown), and the center-line clear- 
ance increased by about 25 percent as shown in figure 9. The initial axial 
disturbance in case 4 was small compared to the angular one. However, the 
amplitude of the axial vibration is large during the first three revolutions as 
can be seen from figure 9. This could be a problem in cases where the ring 
is underdamped and the axial vibration decay slowly. 

Normalized Conir^ 

Figure 10 presents the effect of the normalized coning j3 on the seal be- 
havior, Again from figure 5 we see that case 5, wliich corresponds to 
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/3= 2, is in the unstable regime while cases 3, 6, and 7 are stable. Indeed, in 
figure 10 for case 5 approaches zero while ir ihe other three cases 

Hmin approaches 1.0 (ring parallel to the seat). Comparing the three stable 
cases we see that case 3 with /3= 10 is the most stable while case 7 with the 
relatively large coning /3=100 is the least stable. This result complies with 
the finding of an optimum coning which by equation (39) and for Rj = 0.9 is 
^opt~ Ii'creasing the coning above the optimum not only makes the seal 
less stable but also results in underdamping of the system. Case 7 is severe- 
ly underdamped and as can be seen from figures 10, 8, and 9, the fluctuations 
in cXy, and Z persist much longer than in the other stable cases. 

Dimensionless Group (rj^/rgp)^(Pg - Pi)rQ/K*C^ 

— ■ — " - ft ft 

Figure 11 shows the effect of the group of parameters (r /r_„)'^(p - p.)r^ 
/K*C^. All the four cases compared in figure 11 correspond to a coning /3 = 10 
and, as can be seen from figure 5, with the exception of case 1, are in the 
stable regime of operation. Thus, increasing the pressure differential, for 
example, can make an unstable seal to become stable. As shown in figure 11, 
tliis stabilization effect is associated with some underdamping but both the 
initial disturbance and the fluctuation in die out very rapidly. 

Radius Ratio r^/r^ 

The effect of R^ on stability is clearly presented by figures 4 to 6 from 
which it is seen that wider seals have larger range of stable operation. How- 
ever, from equation (39) we see that for R^ >0.5, which is the practical case, 
^opt decreases as R^. decreases, and while for R. =0.98 is about 102, 
for R. =0. 8 /3^p^ is only 12.5. This has an effect on the critical damping of 
the system. In trying to keep as small as possible to reduce leakage, a 
wider seal may be operating with /3 > /3^p^ and, hence, become underdamped 
and less stable. 
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From the results of the 10 representative cases shown in figures 7 to 11 
we see that, in general, reducing the ring mass m*, its radius of gyration 
r^, and the operating speed o) have a stabilizing effect. Increasing the 
coning angle (3* or reducing the center-line clearance stabilize the 
seal as long as the normalized coning i^+r^/C^ does not exceed the optimum 
coning given by equation (39). If the normalized coning /3 is too high, 

the system is underdamped and tlie range of stable operation is reduced. 
Increasing the pressure p^ on the outer diameter of the seal is beneficial 
from the stability standpoint although with increasing pressure the seal may 
again become underdamped. Increasing the sprii^ constant K* and spring 
location rgp is also favorable for stable operation. An external moment, al- 
though not necessarily affecting the stability, results in a relative tilt between 
the ring and the seat associated with an increase in the center- line clearance 
C and reduction in minimum clearance This, in turn, increases the 

leakage across the sealing dam and can also make the minimum clearance 
dangerously small. A major source of external moment is an uneven spring 
support, hence, effort should be made to eliminate such source. This demand, 
combined with the requirement of stiffer springs for more stabili" r, may call 
for a new approach in designing the flexible support for the seal ring. 

The present analysis treats an ideal seal model in which the seat is per- 
fectly aligned with the shaft, and in which secondary seal friction is neglected. 
In reality, some runout of the seat is almost unavoidable and friction cannot 
be totally eliminated. These two factors, however, do not alter the basic be- 
havior of mechanical face seals, and the three modes of stable, transition, 
and unstable operation discussed in this work were experimentally observed 
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in cases where runout and friction do exist [24]. In reference [22] it is shown 
that the amount of runout does not affect the stability limits of the seal. 

Hence, a stable ring will synchronously track the misaligned rotating seat, 
an unstable ring will contact the seat because of increasing vibration, and in 
the transition region half frequency wobble is superimposed on the synchron- 
ous tracking motion of the rir^. Friction may help in cases of underdamped 
systems but will not allow parallel tracking whenever runout exists. Parallel 
tracking is essential for minimum leakage and from this standpoint friction is 
undesirable and should be minimized. 

CONCLUDING REMARKS 

The dynamic behavior of a noncontacting face seal is analyzed taking into 
account various design parameters and operating conditions. The analysis is 
based on an analytical solution of the Reynolds equation for the fluid film pres- 
sure combined with a numerical solution, using time integration techniques, 
of the equations of motion of the flexibly mounted seal rii^. 

Three modes of seal operation are found that depend on two dimension- 
less groups of parameters. The seal can be either stable, unstable, or oper- 
ate in a transition mode. In the stable regime any disturbance of the ring 
from its equilibrium position decays after awhile. In the unstable regime 
seal failure occurs due to face contact. In the transition region the ring 
wobbles with a constant amplitude at frequency that equals half the shaft 
speed. 

Stability maps for several radius ratios are presented and the effect of 
various parameters on the dynamic behavior is demonstrated and discussed 
through some representative cases. An empirical expression for stability 
threshold of face seals is given along with an expression for optimum coning. 
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It is believed that the present approach will result in a better understanding 
of the mechanism of operation of mechanical face seals, and will enable the 
design of more predictable and reliable seals. 
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APPENDIX 1 

CLOSING FORCE AND RESTORING MOMENT Mgp 

Considering only the support stiffness and neglecting friction due to the 
secondary seal and antirotation locks, the flexible support can be repre- 
sented by a spring constant K*. Assuming a continuous spring distribu- 
tion the spring constant per unit of circumferential length is K*/27rr 
where r is the radius of springs location. Initially the ring is pressed 
against the seat by a spring preload F*_„,.. As the pressure differential 

opi 

is applied and the ring lifts off, the spring force, tending to close the seal- 
ing gap, increases to 


F* = F* . + K*C 
sp spi 


(40) 


The spring force F* in equation (40) is not affected by the nutation y* 

bp 

(see fig. 2) since the additional compression of the springs over one half 
of the ring circumference is counter-balanced by equal relief of compres- 
sion over the other half. 

The total closing force is combined of the spring force F* and the 
outcome of the pressures actir^ on the back side of the ring. If r^ is the 
balance radius (the radius at which the secondary seal contacts the primary 
ring in fig. 1), then the total closing force is 


T* = F|p + u (r 


p„ + ® - rr )p 




Using the seal balance ratio B = ir~ - r, 


fj< 


2 2 

r - r. 
o 1 


= *’lpl + K’C + ^ (r2 - rf) [p^B + p,(l - B) 


the closing force is 
(41) 
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The restoring moment generated by the support is not affected by axial 
displacement but by the nutation 7"* alone. Hence, 


M* = — 
"P 27rr 



27T 


sp'^0 


cos^e d0 = — r^ y* 

&p 2 


( 42 ) 


Normalizing equations (41) and (42) the dimensionless force and moment are 


= *’spi + K ^ (1 + Z) - (l - Ri) [Po® + Pi<l - 


( 43 ) 


M =S r2 

o 


(44) 


In the absence of an external moment the primary ring is balanced at 
Z = y = 0 by the closir^ force F* and the axial force generated by the fluid 
film pressure. From equation (19) it is clear that at equilibrium, Z = y = 0, 
the pressure component = 0. Hence, only the hydrostatic pressure Pg 
contributes to the fluid film axial force F*. At equilibrium when Z = 0 we 
also have 6 = /3 and the fluid film axial force is (ref. 6) 


F* = Trr (r r. ) 


p . + + = — (p„ - p. ) 

2 + ^(1-r/o ' 


(45) 


Normalizir^ equation (45) and substituting for the closii^ force which corre- 
sponds to Z = 0 in equation (43), the closing force F^ becomes 
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F =K~2z+7tRj^(1-R.) 


P. + + 


/3(1 - Rj) 


(P. 


• « * ^ ’ . > rv 

' ° 2 + (3(1 - Rj) 



(46) 
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APPENDIX 2 

THE DESIGNED CLEARANCE C^ AND BALANCE RATIO B 

0 

In analyzing noncontacting seal dynamics the design clearance C^ is 
to be known. It is used in the normalization of various parameters and is 
very important for leakage prediction. 

Using the equilibrium fluid film axial force of equation (45), substituting 
C = Cq in equation (41) and equating (41) and (45), noting that r^ + r. = 2r^^, 
we have at equilibrium 


■p* . 

spi 


2irr(r^ 
m' o 







m - Rj) 
27 i 3(1 - Rj) 



Pl>- 


K*C, 


(47) 


Hence, C^ can be found for any given set of the parameters B, 

(Pq - p^, /3*, K*, r^, and r^,. Alternatively the equilibrium clearance C^ 
can be preselected from leakage consideration, for example, and other pa- 
rameters like the initial spring load F|p., or the balance ratio B, neces- 
sary to maintain the selected clearance can be calculated. 

In noncontacting face seals the range of balance ratio values is limited 

2 

by the normalized coning jS. Dividing equatic^n (47) by Sr^ and rearranging 
we have 



.R„(l 





Ml • Rj) 

2 + Ml- Rj) 


- 2 




(48) 
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For practical reasons the initial spring load has to be positive, hence, 

the right hand side of equation (48) must also be positive. Since the spring 
constant K is always positive, the balance ratio B will have an upper or 
lower limit depending on the magnitude of (P^ - Pj). If (P^ - Pj) > 0 then 
B has an upper limit given by 



/3(1 - Ri) 

2 + /3(1 - R^) 


(49) 


resulting in a balance ratio that is always less than 1.0. If, however, 
(P^ - Pji^) < 0 then B has a lower limit given by 


B 



/3(1 - Rj) 

2 + /3(1 - Rj) 


(50) 


which results in a balance ratio that is always larger than 0. 5. Equation (48) 
with the requirement of F . > 0 along with the conditions (49) and (50) 

bpi 

provide the range of balance ratio and spring constant values from which K* 
and B can be selected for a desired equilibrium clearance C^. After K* 
and B are selected to meet these conditions, the required initial spring 
load can be calculated from (47). 

It is worthwhile noticing that in noncontacting face seals, the required 
balance ratio B, as given by equations (49) and (50), may be negative when 
p^ > p^, and may be larger than unity when p^ < p^. The first case corre- 
sponds to r^j > r^ while the second corresponds to r|^ < r^. 
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TABLE I. - CHARACTERISTICS OF REPRESENTATIVE 
TEST CASES, R. = 0.9 


Case 

number 

(rg/rgp)^m*a)‘VK* 


X 

M /Sr? 
e' o 

1 

32 

10 

90.25 

0 

2 

21.95 

10 

90.25 

0 

3 

8 

10 

90.25 

0 

4 

8 

1 

10 

90,25 

0.01 

5 

8 

2 

90.25 

0 

6 

8 

50 

90.25 

0 

7 

8 

100 

90.25 

0 

8 

32 

10 

361 

0 

9 

32 

10 

902.5 

0 

10 

32 

10 

3610 

0 
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Figure 8. - Variation of normalized tilt with time for seals operating at various 
stability modes, R; ■ 0.9. 
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Figure 9. - Variation of dimensionless axial displacement with time for seals op- 
erating at various stability modes, R| * CL 9. 




